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ABSTRACT 

It has recently been noted that many discs around T Tauri stars appear to com- 
prise only a few Jupiter-masses of gas and dust. Using millimetre surveys of discs 
within six local star-formation regions, we confirm this result, and find that only a 
few percent of young stars have enough circumstellar material to build gas giant plan- 
ets, in standard core accretion models. Since the frequency of observed exo-planets 
is greater than this, there is a 'missing mass' problem. As alternatives to simply ad- 
justing the conversion of dust-flux to disc mass, we investigate three other classes of 
solution. Migration of planets could hypothetically sweep up the disc mass reservoir 
more efficiently, but trends in multi-planet systems do not support such a model, and 
theoretical models suggest that the gas accretion timescale is too short for migration 
to sweep the disc. Enhanced inner-disc mass reservoirs are possible, agreeing with pre- 
dictions of disc evolution through self-gravity, but not adding to millimetre dust-flux 
as the inner disc is optically thick. Finally, the incidence of massive discs is shown to 
be higher at the protostellar stages, Classes and I, where discs substantial enough to 
form planets via core accretion are abundant enough to match the frequency of exo- 
planets. Gravitational instability may also operate in the Class epoch, where half the 
objects have potentially unstable discs of >30 % of the stellar mass. However, recent 
calculations indicate that forming gas giants inside 50 AU by instability is unlikely, 
even in such massive discs. Overall, the results presented suggest that the canoni- 
cally 'proto-planetary' discs of Class II T Tauri stars have globally low masses in 
dust observable at millimetre wavelengths, and conversion to larger bod- 
ies (anywhere from small rocks up to planetary cores) must already have 
occurred. 

Key words: circumstellar matter - planetary systems: protoplanetary discs - plan- 
etary systems: formation 



1 INTRODUCTION 

It is often assumed that the discs of classical T Tauri stars 
are those in which planets form. These 'Class IF objects 
have disc-dominated dust emission in the infrared to mil- 
limetre, in contrast to fainter remnant discs around Class III 
(weak-line) T Tauri stars, or the envelope-dominated emis- 
sion of the earlier-stage Class and I protostars. Provided 
that dust persists for a few Myr in Class II discs, along 
with the gas com ponent that contains the bulk of the mass 
l|Dent et al.ll2005l ). then gas giants should be able to form 
by core accretion followed by accumulation of a t hick atmo- 
sphere l|Pollack et al.lll99d : fHubickvi et al.ll2005l ). However, 
it has been noted recently that the masses of T Tauri discs 
as derived from millimetre observations are typically small 
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compared to the reservoirs n eeded for gas-giant formation. 
And rows fc Williams! i|2007af ) inferred that the median Class 
II disc comprises only 5 Jupiter masses of dust and gas, from 
deep surveys of the Taurus and Ophiuchus regions. In con- 
trast, at least 20 Jupiter masses was present in the Minimum 
Mass Solar Nebula (MMSN), on the b asis of the r ocky con- 
tent now locked up in planetary cores (|Davis!l2005l ). In mod- 
els, a mass-reservoir a few times larger is actually needed for 
core-accretion to successfully build the the Sun's giant plan- 
ets. For example. iDodson- Robinson et all (|2009l ) find that a 
120 Jupiter-mas£] or 6-MMSN disc can make all four giants. 
There are also examples of more massive exo-planets than 
the Sun's giants, probably requiring more substantial discs 



1 The convenient relation Mj„ piter ) f» Mg un /1000 is adopted 
throughout, rather than the true inverse mass ratio of 1047. 
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from which to form. Hence, if exo-discs only comprise a few 
Jupiter-masses, there is a 'missing mass' problem at these 
early stages. 

Further, the surface densities that are the critical pa- 
rameter in core accretion models are typically measured to 
be lower than expected, for example in relation to form - 
ing Jupiter at around 5 AU j Andrews fc Williams! l2007rj ). 
Millimetre dust emission is a good choice for measuring to- 
tal disc-masses, being of low optical depth and reasonably 
bright, but limited angular resolution means that material at 
a few AU is not resolved. Although the bulk of the disc mass 
should be at large scales l|Wvatt et al.ll2007l ). discontinuous 
profiles with a strong central mass peak can not be ruled 
out. However, if th e profiles are smooth, the ty pical outer 
radii of « 200 AU (| Andrews fc Williams! l2~007bD are much 
larger than the Solar System, exacerbating the missing-mass 
problem by further reducing the inferred inner-disc surface 
densities. 

Several types of solution to the missing-mass problem 
have been proposed. The conversion of millimetre flux to 
dust mass is via the grain temperature and emissivity, and 
while the former can be estimated from the wavelength of 
peak emission, the latter is m ore uncertain and comes from 
modelling of grain properties (|Drainej|2006l ). Decreasing the 
adopted emissivity would enhance the derived disc masses, 
and this could be reasonable for a different composition, 
or a size distribution where much of the mass is in large 
grains. Such 'pebbles' are not conspicuous in the millime- 
tre meas urements but can be d etected by radio continuum 
emission (jRodmann et al.ll2006T l. In addition, the total disc 
masses could be wrong if the assumed gas-to-dust mass ra- 
tio of 100, i.e. inherited from the interstellar medium, is in 
error. A higher conversion factor would be correct if some 
of the dust has evolved into non-observable forms (any- 
where from small pebbles up to planetesimals) but 
the gas is still all in the disc. The true disc masses should 
then be higher if a more appropriate gas-to-dust ratio were 
used. Similarly, if some of the gas has also been accreted 
into planets, the disc mass estimates are even more un- 
certa in but will underestimate the init ial reservoir. Gener- 
ally, lAndrews 8z Williams! (|2007ah have shown that 
the disc masses estimated via the amount of accre- 
tion onto the star are substantially higher than the 
masses derived directly from the observed dust, con- 
firming that there is indeed a discrepancy within the 
available data. 

There are however some difficulties with these kinds 
of solutions to the missing-mass problem. Firstly, a bulk 
increase in the adopted emissivity w ill affect all systems, 
and t he few discs of > W0Mj up u er l| Andrews fc Williams! 
2005, 2007a. b) would then be comparably massive to the 
host star. This situation is likely to be dynamically unstable, 
i.e. such a disc would tend to fragment to form a compan- 
ion star or brown dwarf, and not remain to be observed in 
the Class II era. Secondly, where larger grains are detected 
in long-wavele ngth data, the d e rived disc masses do not al- 
ways increase. iRodmann et al.l (|2006h found systems where 
the mill imetre spectr al index is small - a signature of grain 
growth |Drainel 2006) - with three discs having indices as low 
as 0.5-0.8 in the 1-7 mm regime. However, masses derived 
from their 7 mm fluxes range from 20 % lower to 70 % higher 
than the estimates based on 0.85 mm data, rather than being 



substantially increased. Thirdly, the bulk H2 gas component 
has in some cases been detected by infrared spectroscopy, al- 
though this is only sensitive to warm layers and so the total 
gas reservo i r mu st again be estimated. In their models 
iBarv et al.l (|2008T ) find that only a few Jupiter masses of hy- 
drogen at ~ 10 — 30 AU is needed to reproduce line fluxes, 
but many systems a few Myr old are not detected, so pre- 
dominantly high gas-to-dust mass ratios are not supported. 
The gas could already have bee n accreted into plan etary 
atmsopheres in favourable cases (|Hubickvi et al. I l2005h . but 
it may instead have been simply lost by m echanisms such as 
photo-evaporation (Erc olano et alj[2009h . 

Here we examine some alternative solutions to the 
missing-mass problem. We first compile a millimetre-based 
distribution of Mdisc over several regions of young stars, 
and consider if there are globally enough high-mass discs to 
match the fraction of Sun-like stars observed to host gas gi- 
ants. Finding that there is still a disc-mass deficit, we then 
examine three alternative classes of solution: that the mass 
reservoir could be used more efficiently if planets migrate 
while forming; that mass is hidden at small radii because 
it is optically thick; and that the mass has disappeared be- 
cause it is accreted into planets very early on, even at the 
protostellar stage. 



2 DISC MATERIAL 

Star formation near the Sun is well characterised by the 
molecular clouds of the Gould Belt, lying at distan ces of 
around 150 to 500 pc (|Ward-Thompson et al.l 120071 1. The 
belt includes at least a dozen regions of young stars, of which 
six have deep millimetre-wavelength surveys for dust discs. 
Other Gould Belt regions have few millimetre-detected discs, 
either because the stars are still mostly very deeply embed- 
ded so that intensive millimetre interfometry is needed to 
identify compact structures (e.g. Serpens, Corona Australis, 
IC 5146), or the large cloud areas have yet to be thoroughly 
mapped (e.g. Cepheus Flare, Pipe Nebula). The disc sur- 
veys included here are listed in Table 1, and cover the Tau- 
rus/Auriga, p Ophiuchus, Chamaeleon and Lupus regions, 
the Orion Nebula Cluster and the IC 348 cl uster in Perseus. 
The median stellar ages in these clouds (|Palla fc Stahlej 
2000) are 1-2 Myr, although the disc surveys include Class O 
and I protostars that may be 1 Myr old. The surveys have 
target selection biases, and in particular faint Class III sys- 
tems may be under-represented. Also, not all the young stars 
were observed, with surveys sampling ~ 30— 150 targets out 
of ~ 150—1500+ objects in each cluster. However, these data 
are the best snapshot presently available of the properties 
of potential planet-forming discs. The environments range 
from sparse star formation in Tau/Aur to the dense ONC, 
with the total star-count greatly dominated by this latter 
region; hence, its discs are of particular interest as they may 
represent archetypal planet-forming conditions. 

Subsequent discussion implicitly assumes that these 
discs we observe today are equivalent to those that formed 
planets in the past - typically a few Gyr ago when consider- 
ing the Sun-like hosts of observed exo-planets. We also as- 
sume that all the star-formation environments that existed 
in the past are also observed today; that the exo-planets ob- 
served are not from a peculiar subset of stars; and that the 
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Tabl e 1. Summary of disc masses from millime tre surveys. Cluster star counts are from lPorras et al except for Lup l|Merfn et al] 

2008) and Tau/Aur jKenvon fc Hartmannlll995tl , and may not be identical regions to those surveyed for discs. Ages given are the medians 
within each cluster from lPalla fc Stahlerl l|2000l l. Data at A mm are from interferometers for ONC and IC 348 and otherwise from single 
dishes — these photometry results may include cloud emission in crowded regions. The number of discs detected is given as a subset of 
the total number of targets including all spectral types and evolutionary stages; in IC 348 a net detection only was made over all the 
stars observed. Subsequent columns list the fraction of MMSN discs (above 20 Mj„ p ) and Poisson errors, the median mass (approximate 
where over half the data are upper limits), and the 100 %-complete mass limit of the survey. 



cluster 


age 


N stars 


A 


N observed 


Immsn 




(Myr) 




(mm) 


(Ndiscs) 


(%) 


ONC 


1 


1471 


0.88 


55 (26) 


11 ± 4 


Oph 


1 


337 


0.35-1.3 


147 (64) 


18 ± 4 


Tau/Aur 


1.5 


289 


0.35-1.3 


153 (94) 


17 ± 3 


IC 348 


1.5 


345 


3 


95 (0) 


~ 


Lup 


2 


159 


1.3 


32 (12) 


6 ± 4 


Cha 


2 


180 


1.3 


36 (16) 


8 ± 5 



mode of planet formation that dominated in the past has 
not greatly altered as for example the Galaxy has become 
more enriched in heavy elements. 

2.1 Millimetre-based disc masses 

Table 1 summarises the results of the these millimetre disc 
surveys, and presents global properties for all the systems 
observed in each region. The results are adopted directly 
from the survey papers without attempting to align slightly 
different assumptions about dust temperatures and emissiv- 
ities. All the surveys were complete down to an MMSN of 
20 Mj up , apart from that of IC 348 where this represents a 
2.5<t upper limit for each star. 

With this very comprehensive dataset of 518 stars, it 
is clear that the fraction hosting an MMSN when observed 
is small. No cluster has an incidence /mmsn > 20 % and 
the overall average is /a/a/sjv = 12 ± 2 % (Poisson er- 
ror). The median disc of each cluster comprises only a few 
Jupiter masses, down to as little as 1 — 2Mj up in the deepest 
data. There are no marked differences between sparse and 
dense clusters, and }mmsn does not decline significantly 
with time. All of the clusters are within ~ 1.5<r of the global 
/mmsn, except for IC 348 where no MMSN are seen (but 
the survey is barely sensitive to this depth). 

Figure 1 shows the estimated underlying distribution of 
disc masses. The discs have been divided among bins cov- 
ering a factor of 3 in mass, with the highest three bins rep- 
resenting an MMSN or more. The raw counts were scaled 
up to infer the numbers of discs that would be expected 
for all the stars known in the cluster (Table 1), and these 
counts are then stacked to show the estimated contribution 
of each cluster. For uniformity the upper limits were treated 
as values, so the plotted distribution is the most optimistic 
one from the view of planet-forming potential. Two clusters 
were treated as s pecial cases: the ONC, where complete- 
ness information dMann fc Williams! [2009) was used to in- 
fer limits, and I C 348, where the nominal measured fluxes 
(|Carpenterll2002l ) were treated as data points where positive, 
and the mean 2 Mj up disc was adopted for all other starfl 

2 This avoids the problem of logarithms of negative values but is 
not strictly flux-conserving, with the mean signal per disc being 
inflated to the +4<r bound. 
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Figure 1. Distribution of disc masses from surveys listed in Ta- 
bic 1, scaled by the number of stars in each cluster and treating 
upper limits as values, hence maximising the Mdisc bin they lie 
in (see text). The three right- most bins correspond to one MMSN 
upwards. On the left side, sparse data for low-mass discs (7 % of 
stars) are omitted for clarity. 

For Tau/Aur and Oph only, some limits or detections below 
2 Mj up were reached, but these are omitted from the plot 
for clarity. 

Figure 1 also implies that the number of discs with 
MMSN-like planet forming potential is low. Including the 
scalings by total number of stars per cluster, /mmsn is 
now 11 %, very similar to the uncorrected value above of 
12 %. The proportion of discs in the top two bins, above 
the thresholds in many gas giant formation models, is only 
about 0.5 %, and such massive discs were only found in the 
Tau/Aur and Oph regions. 

2.2 Disc population 

Survival analysis is used next to make a better estimate of 
the underlying distribution of disc masses, given the vary- 
ing limits between and within surveys. The Kaplan-Meier 
estimator is appropriate when the limit obtained is inde- 
pendent of the mass of the disc being observed, a condition 
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log Mdisc (Msun) log Mdisc (Msun) 

Figure 2. Kaplan-Meier estimates of the underlying disc population (see text). The upper and lower stacked histogram bars indicate 
limits and detections in each mass bin respectively, and the line is the K-M fit. The right-most three bins comprise 1 MMSN upwards. 
Left: weighting each cluster equally; Right: weighting each cluster proportional to its total number of stars from Table 1. 



that is reasonably well met here. The K-M tests were per- 
form ed on the log-disc-m ass data using the ASURV pack- 
age (jLavallev et al.lll992l ) and the results are shown in Fig- 
ure 2. In the first panel, all of the 6 clusters of Table 1 were 
weighted equally, by taking 32 stars from each (i.e. all of 
those observed in Lupus and a random set of 32 from each 
of the other clusters) . In the second panel, the clusters were 
weighted according to the total number of stars in each, i.e. 
taking all 55 observed discs from the ONC and a proportion- 
ally smaller random set from the other clusters (down to 6 
stars in Lupus, i.e. 55 x {N Lup = 159)/(N nc = 1471)). 
These two approaches attempt to take into account a dif- 
ferent mix of star-forming environments that could have oc- 
cured at past times. 

In both cases, Figure 2 shows a broad distribution of 
masses plus a large tail of very insubstantial discs. The 
main distribution peaks at 5 Jupiter masses, while about 
30 % of stars have such negligible discs (left-most bin) that 
the dust mass is < IMq - these might not even be able 
to make a single terrestrial planet. At the high-mass end, 
weighting all the clusters equally yields an MMSN fraction 
of 13.5 %, while the -/V star -weighted calculation gives 10 % 
MMSN. These results are very similar to the previous re- 
sults of Jmmsn ~ 11 — 12 %. Hence, even if the modes of 
star formation differed in the past, the fraction of discs ca- 
pable of forming gas gaint planets appears to have been low, 
and many stars seem to have been incapable of any planet 
formation. 

2.3 Disc mass evolution 

It has already been suggested that higher disc masses at 
earlier stages mi ght be more conducive to plan et formation, 
and for example Andrews fc Williams! l|2007al ) have shown 
that this is the case for Taurus and Ophiuchus. There the 
median disc mass is 5 Mj up for Class II stars but 15 Mj up 
for the precursor Class I objects. Table 2 gives a similar 
summary of the data for all six star-forming regions of Ta- 
ble 1. No weighting of the clusters is made, as wide varia- 
tions in detection-rates and proportions of different classes 



Table 2. Disc masses sub-divided by object class. Data from all 
6 regions of Table 1 are added together without weighting. There 
are 38, 25 and 339 objects in Classes I, FS and II respectively, and 
the corresponding fractions of upper-limit data are 5, 24 and 43 %. 
For each pair of lines, the upper includes detections only, and the 
lower (bracketed) includes the upper-limit systems but treated as 
zero disc-mass. The Poisson error is given for the frequency of 
MMSN. 
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Flat Spectrum 


Class II 


mean M disc 


41 


14 
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(M Jup ) 


(39) 


(11) 


(8) 


median M^ sc 


« 15 


10 


9 


(M Jup ) 


(« 15) 


(10) 


(2) 


}mmsn 


47 ± 11 


26 ± 12 


22 ± 3 


(96) 


(45 ± 11) 


(20 ± 9) 


(13 ± 2) 



observed make this difficult to do systematically. Where ob- 
jects are not explicitly classifi ed, we treat as all-Class I I both 
the ONC proplyds (following iMann fc Williamsl l|2009l )) and 
the slightly older IC 348 stars (for consistency, although 
in fac t the proportion o f II:III objects here is about 
30:70 |Muench et alj|2007l )). Class III stars are not included 
in Table 2 because there are very few dust detections and 
so average properties are not meaningful. 

Table 2 confirms that the fraction of MMSN is higher 
for Class I objects than at either of the subsequent Flat- 
Spectrum or Class II stages. The median disc mass is below 
an MMSN at all stages, but at Class I only just sc0, and also 
7/38 objects (18 ± 7 %) have disc-masses of 2-25 MMSN as 
in many giant-planet formation models. The previous re- 
sults on the low fraction of MMSN are clearly dominated by 
the preponderence of Class II discs, with Class I/FS objects 
representing briefer stages and thus fewer examples seen at 
snapshots in time. The mean disc-mass declines after Class I, 
but the median mass only drops after the FS stage, when a 

3 Millimetre photometry may somewhat over-estimate Class I 
disc masses, if the remnant envelopes contribute some flux. 
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large tail of upper limits appears within Class II. Combined 
with the remnant discs of Class III, many systems have very 
low masses by the few-Myr epoch. 



2.4 Comparison to planetary systems 

We now consider the mass budget needed in early-time discs 
to form t he observed extrasolar planets. Analysis of Doppler 
surveys (jCumming et alJ l2008h shows that approximately 
18 % of Sun- like stars should host a giant planet, of above 
about Saturn's mass and orbiting within 20 AU. The re- 
sults are based on extrapolations of detections made so far, 
but are consistent with different extrapolation methods. The 
typical planet mass is about that of Jupiter, and a little 
over half of these systems should host a body orbiting at 
about 3 AU outwards, comparable to the location of the 
Sun's gas giants. Within this region, microlensing surveys 
indicate that Neptune-mass pla nets are at least t hree times 
more abundant than Jupiters (|Sumi et alj |2010). Multiple 
planetary systems are expected to be common, with over a 
quarter of those found by Do ppler wobble alread y known to 
have two or more gas giants (|Wright et al.ll2009l ). 

If giant planets orbiting at a few AU are common, early- 
time disc-mass reservoirs of a few MMSN are expected. A 
high surface density of solids is needed to build massive plan- 
etary cores, with subsequent accr etion of large masses of 
volati les in the case of gas giants. iDodson-Robinson et al.l 
l|2009f l find that with 120 M Jup of material (6 MMSN), all 
four of the Sun's giant planets could be made, including 
non-water ic es for speeding up the sticking of planet cores. 
iDeschl (|2007l ) argues for a steeper-profiled initial nebula with 
later outwards migration of the ice giants, but still finds at 
least of order 60 Mj uv (3 MMSN) within 30 AU is required^. 
lAlibert et all (|2005h have accounted for various properties of 
Jupiter and Saturn in model discs of 35-50 Mj up , or 1.75-2.5 
MMSN. 

According to thse models, giant-planet systems should 
be rare, if the discs seen in the time-snapshot from the 
millimetre surveys represent true budgets of planet-forming 
material. Only 4 % of stars host discs of > 1.75 MMSN, 
as required in the most optimistic models for gas giants 
l|Alibert et al.l [2005) . However, the giant planet detection 
rate is already 7 % in the regime where surveys are close 
to complete, and ext rapolated to 18 % out to 20 AU 
i|Cumming et al.ll2008l), or even high er including lower-mass 
'ice giant' planets 1 Sumi et alJ 2010). Therefore, even if only 
the most massive discs go on to form planets, there is ap- 
parently a major discrepancy in the statistics, with more 
exo-planet systems than discs suitable to make them. 



3 SOLUTIONS 

As discussed above, the millimetre dust data may under- 
estimate the true total disc masses, if any of the several 
assumptions are incorrect. Corrections to dust emissivities, 



4 So me of these authors' as s umpt ions are slightly different to 
ours . IDodson-Robinson et al,l adopt a gas-to-dust mass ra- 

tio of 70 not 100. and IDeschl <2007ft adopts a smaller MMSN of 
13 Mj up for a steep disc profile. 



relative populations of large and small grains, and gas-to- 
dust proportions could globally shift the distributions to- 
wards higher ma sses. However, to repro duce « 18 % of giant- 
planet systems (ICumming et al.l I20 08H requires all discs in 
the millimetre surveys of nominally > lOMj^p to be planet- 
formin g, whereas the mod els need at least 35 Mj up of ma- 
terial (lAlibert et al.ll2005h . Since an increase of a factor of 
~ 3.5 to all millimetre-derived masses would be both sub- 
stantial and ad hoc, we now consider if other assumptions 
may be wrong. In particular, we consider further solutions 
to the missing-mass problem, of three kinds: 

(i) mass is used more efficiently than was the case in the 
Solar System, with migrating planets sweeping up much of 
the disc; 

(ii) more mass is 'hidden' at small scales in the centre of 
the disc, with no increase of millimetre emission as the dust 
is optically thick; 

(iii) the T Tauri discs observed are only remnants, with 
most of the dust grains growing into larger bodies during 
the protostellar stage — this could range from cm-sized 
rocks up to fully-formed planets under the disc in- 
stability hypothesis, discussed below. 



3.1 Migration 

The Solar System gas giants are thought to have 
migrated only slightly f rom their formation positions 
i|Minton fc Malhotra|[2009l. e.g.), but exo-planets may have 
undergone larger orbit changes. 'Hot Jupiters' must have mi- 
grated significantly inwards, potentially sweeping large disc 
regions clean of material. Highly efficient usage of the disc 
mass-reservoir - in co ntrast to < 10 % of the MMSN that 
ended up in planets (|Davis! 120051 ) - might hypothetically 
explain how the discs can be of low-mass and yet still form 
planets. However, long formation times in the outer disc will 
set limits, as cores here would accrete mass slowly and have 
migrated little by the time the disc disperses. Therefore it 
is implausible that the entire disc could be swept up into 
planets. 

If disc-sweeping during migration is important, there 
should be observational signatures in the properties of exo- 
planets. For example, in systems with two planets, the 
outer one would tend to be of lower mass if it has passed 
through areas already swept clean by the migration of the 
inner body. If the two planets have both ended up close- 
in to the star, this trend should be particularly marked, 
whereas if the two planets are well separated, they could 
have swept the inner and outer disc independently and 
so not show such a trend. We searched for such signa- 
tures among multiple-planet systems, using only the pa- 
rameter space Mpianetsini 0.3A/j up and semi-major axis 
a ^ 3 AU, where discoverie s are thought to be almost com- 
plete l|Cumming et al 1l2008l ). This avoids a bias where lower- 
mass planets on large orbits are harder to detect (i.e. smaller 
Doppler-wobble) and so such systems would be missing from 
the multi- planet plot. Considering the well-studied volume 
of Sun-like (FGK) stars within 60 pc, 10 systems have two 
planets within our chosen M, a parameter space and v And 
has three (Figure 3). 

These data do not support the idea that accreting plan- 
ets have swept up much of the disc. In particular, only one 
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Figure 3. Observed masses and semi-major axes for systems with 
two or three planets (from http://exoplanct.cu), in the regime 
where searches are close to complete (see text). Lines connect 
planets within each system, as named at the right. 



system in Figure 3 (HIP 14810) appears to have a 'starved' 
outer planet and both bodies ending up near the star. Also, 
5/11 of the systems have outer planets that are more massive 
than the inner planet, contradictory to the idea of sweeping 
through already-depleted regions of the disc. 

From a theoretical perspective, a further problem with 
the idea that migration can aid planet-building is that most 
of the gas accretion occurs in a runaway phase, with a 
timescale much shorter than typical migration timescale s 
l|Pollack et al.lll996l ; llk"oma et aljlioool ; iBrvden et al.ll2000l ). 
Migration is thus implausible as a mechanism that allows 
a planet to sweep up its bulk in gas from a large region 
of the disc. However, m igration could aid in the growth of 
the s olid planetary core (jHourigan fc W ard 1984; [Rice et al.l 
2003), and this could allow the core to reach a critical mass 
to attract an atmosphere, while there is still sufficient gas 
in the disc. While migration is thus still important to plan- 
etary evolution, neither theory or observational constraints 
suggest that it solves the 'missing-mass' problem of discs. 



3.2 Centrally-concentrated discs 

Small corrections for opacity are known to be needed when 
convert ing millimetre dust emission f rom T Tauri discs into 
masses. lAndrews fc Williams! l|2007bl ) estimate that the ra- 
tio of optically thick to optically thin submillimetre emis- 
sion is typically around 0.3, as more opacity would result 
in flattening of the spectrum. However, if the discs have a 
massive central inner region, on unresolved scales of tens of 
AU or less, then this could be much more optically thick 
- cont ributing significan t mass, but little e xtra m illimetre 
signal. IZhu et all l|2009l ); iRice fc Armitagel l|2009h have re- 
cently shown that if discs are massive with respect to the 
star, then transport of angular momentum through disc self- 
gravity does in fact lead to a pile-up of material at smaller 
radii. A quasi-steady-state is reached in which ~ 80 % of 
the disc mass ends up within 10-20 AU of the star, with a 
drop to lower surface densities in the model outer disc ex- 
tending to 50 AU. There is thus a physical basis for the idea 
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Figure 4. Spectral energy distributions for (top) a 0.05 Mq (50 
Jupiter-mass) disc with a power-law mass distribution, and (bot- 
tom) the same disc but with more mass added inside 10 AU for a 
total of 100 Jupiter masses. The curves are for different viewing 
angles, with the stellar emission strongly blocked by discs that 
are closer to edge-on. 



of a central mass concentration, at scales relevant to planet 
formation. 

Figure 5 shows representative spectral energy distribu- 
tions (SED) of the star plus disc system - generated using 
the H Q-CHUNK: 3D radiation transfer code i| Whitney et all 
2003) - for a standard power-law distribution of disc mass 
(top panel), and a case where more material has been artif- 
ically added inside 10 AU to double the total disc material 
(bottom panel). There is very little difference in the SED of 
the power-law and centrally-enhanced discs, especially in the 
millimetre regime that is canonically 'mass tracing'. In fact, 
the centrally condensed disc has a submillimetre flux that 
is slightly lower than that from the lower-mass, power-law 
disc. Therefore, if real discs are in fact mas sive, a central pile- 
up w ould be both theoretically predicted jRice fe Armitagel 
120091 ) and not detected in millimetre images where the inner 
disc is unresolved. This provides a potential solution to the 
missing mass, at least for discs with moderately high mass 
estimates already. 



3.3 Disc masses at early times 

As discussed above, discs have higher mass at earlier evolu- 
tionary stages. At Class I, the fraction of objects hosting 2- 
25 MMSN is around 18 %, in the regime of Jupiter-building 
models, and also matching the 18 % of gas-giant planet sys - 
tems inferred around mature stars l|Cumming et all [2008 V 
However, if planetary cores can grow this early on - which 
could indeed occur if migration were to accelerate the core 
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Table 3. Properties of Class protostars observed with millimetre interferometry. Disc masses are lowest- or single-temperature results, 
or from models where a range is given (see papers cited). The results shown are for the longest published wavelength, with > noting that 
emission around 1 mm may be optically thick. Bolometric luminosities are values generally adopted in the literature. The stellar masses 
are derived assuming the protostar's luminosity is powered by infall (see text). 



protostar wavelength Mdisc Lbol M*(Lbol) Mdisc/M* reference 

(mm) (M Jup ) (L Q ) (M ) 



L 1448 IRS2 


2.7 


~0 


5 


~0.08 


~0 


Kwon et al. (2009) 


NGC 1333 IRAS4A 


2.7 


~0 


14-16 


~0.24 


~0 


Loonev et al. ("20031 


B 335 


1.2 


(>)4 


3 


~0.05 


>o.i 


Harvev et al. f 20031 


VLA 1623 


0.8 


(>)4-10 


1 


~0.02 


^0.4 


Brown et al. (2000) 


Ser SMM3 


0.8 


(>)4-16 


8 


~0.13 


>o.i 


Brown et al. f 20001 


L 1448 IRS3B 


2.7 


~30 


7 


~0.11 


~0.3 


Loonev et al. f2003) 


L 1551 IRS5 S 


7 


~30 


~7 


~0.11 


~0.3 


Rodriguez et al. ( 19981 


L 483 mm 


2.6 


$40 


9 


~0.14 


$0.3 


Jorgcnsen et al. (20041 


L 1448 C 


1.4 


(>)40 


5 


~0.08 


£0.5 


Schoier et al. (20041 


Ser SMM1 


0.8 


(>)50-80 


16 


~0.7 


>o.i 


Brown et al. (20001 


L 1551 IRS5 N 


7 


~60 


~23 


~0.4 


~0.2 


Rodriguez et al. ( 1998) 


NGC 1333 IRAS2A 


2.7 


~60 


16-30 


~0.4 


~0.2 


Loonev et al. (2003) 


NGC 1333 SVS13B 


2.7 


~60 


~22 


~0.35 


~0.2 


Loonev et al. (20031 


NGC 1333 SVS13A 


2.7 


~80 


~19 


~0.3 


~0.3 


Loonev et al. (20031 


L 1157 


2.7 


~110 


11 


~0.18 


~0.6 


Kwon et al. (2009) 


NGC 1333 IRAS4B 


2.7 


~120 


5 


~0.08 


~1.5 


Loonev et al. (20031 


IRAS 16293 A 


1.4 


(S) 2 40 


~11 


~0.18 


>l-4 


Schoier et al. (2004) 


IRAS 16293 B 


7 


~300-400 


16 


~0.26 


~1.4 


Rodriguez et al. (20051 


Ser FIRS1 


1.1 


(>)1000 


16 


~0.75 


>l-3 


Enoch et al. (20091 



growth as noted above - then in principle the process could 
even start in the Class protostar phase, within the first 
0.1 Myr. At this time the central object is only partly ac- 
creted, and much of the mass of the system is still in a 
circumstellar envelope. Millimetre photometry will be dom- 
inated by the envelope flux, but many of the Class ob- 
jects have now been studied with millimetre-wavelength in- 
terferometry. This resolves out the large-scale envelope emis- 
sion, leaving in most cases a compact or point-like signal 
attributed to a very young disc. 



Table 3 lists these interferometry-based results, in or- 
der of increasing disc mass. The mass estimates were made 
with different methods and assumed temperatures and emis- 
sivities, and may be lower limits at w 1 mm if the inner 
disc has non-negligible opacity, so only rough comparisons 
may be made. However, among these 19 Class objects, 
half now host upwards of 2 MMSN, increased three-fold 
over the equivalent population within Class I. This contin- 
ues the trend of finding more substantial discs at each ear- 
lier phase studied, and is very promising for forming many 
future planetary systems. If these larger dust masses 
at the Class stage represent most closely a true 
'initial mass reservoir' for planet formation, then 
the small quantities observed during Classes I, FS 
and II suggest that dust has been converted into 
larger pebbles and boulders, and perhaps planetes- 
imals if core accretion has p roceded successfully. In 
the g rain growth models of Dull emond &; Dominikl 
(2005), metre-sized boulders are abundant by ages 
of a few 10 years, i.e. in the Class I stage, and plan- 
etesimals can agg regate into planetary cores over a 
similar timescale (Hu bickvi et al] |2005l. 



3.3.1 Planet formation by gravitational instability 

In addition, as the protostar is largely assembled in the 
Class phase, it must for part of this period have much 
less than its final mass. Therefore, as the discs are mostly 
substantial, Mdisc/M sta r can be high, raising the possibil- 
ity of forming planets by the alternative mechanism where 
gravitational instability leads to disc fragmentation. To at- 
tempt to quantify this, we estimate masses for the protostars 
assuming that the bolometric luminosity is supplied by the 
energy released in infall. The accretion rate is taken to be 
10~ 5 Mo yr~\ i,e, building most of an ~ 1 Mq star in about 
0.1 Myr (|Evans et al.ll2009l l. The mass of the protostar M star 

is then derived from L bot (GM s ta r dM/dt) / Rstar, taking 

Rstar to be w 5R Q after lAndre fc Montmerld' l|l994 l. The 
protostellar radii have not been measured but for pre-main- 
sequence phases shou ld be moderately independent of mass 
I Whitney et alj |2004| ). Mdisc/M s t ar values may be uncer- 
tain by factors up to a few, as disc opacities, stellar radii 
and accretion rates are not known for individual systems. 

The estimates in Table 3 suggest that half of the discs 
should be unstable to fragmentation, with Mdi S c/M s tar 
of roughly 0.3 or more. In four cases, the disc ap- 
pe ars to actually out-wei g h the star, as in some models 
of iHueso fc Guillotl |2005l ); iRice et all (I2009T ). The stellar 
masses are low, as expected if the stars are not yet fully 
built. The protostars estimated to be the most massive also 
tend to have higher-mass discs, as do the objects in binary 
systems (i.e. names ending in A, B, N or S). However, most 
of the short-wavelength data happens to be for the single 
systems, so these could simply have disc masses that are 
under-estimated due to dust opacity. 

One problem with the idea of forming planets early 
on via disc fragmentati on is that theoretical calcula- 
tions, both semi-analytic (|Matzner fc Levin 2005; Rafikovl 
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120051; iRice fe Armitage! 120091; IClarke! 120091) and numerical 
l|Bolev et al.ll2006l ; IStamatellos fc Whitworthll20oi ). suggest 
that planetary mass bodies will not form inside ~ 50 AU. It 
is likely that these extremely m assive Class discs wo uld in- 
stead become globally unstable l|Lodato fc Ricell2005h . mov- 
ing large amounts of mass to large radii, where fragmenta- 
tion c ould produce substellar, or even stellar, mass compan- 
ions |Stamatellos fc Whitworthl |2009h ■ However, the pres- 
ence of spiral structures in the massive disc could help to 
collect together solid particles (lHaghighipour fc Boss|[20o3 ; 
IRice et all \2004 ) and thus accelerate planet formation via 
the mechanism of core accretion. Thus massive discs at very 
early times may still be a promising signpost to abundant 
planetary systems around mature main-sequence stars. 



4 CONCLUSIONS 

The compilation of millimetre data for circumstellar discs 
in 1-2 Myr-old star formation regions confirms that the ob- 
servable mass reservoirs are small, at these snapshots in 
time. The fraction of discs apparently capable of form- 
ing gas giants by core accretion is generically less than the 
fraction of observed exo-planet systems. An early start to 
planetary growth is favoured, as the fractions of exo-planets 
and suitable proto-planetary discs do match in the Class I 
stage. We also considered other solutions to the 'missing 
mass' problem. Sweeping up the disc material efficiently into 
planets has both observational and theoretical difficulties. 
More promisingly, viscous evolution through disc self-gravity 
could produce a central concentration of mass that would be 
optically thick even in the millimetre, and so essentially in- 
visible. At very early times (Class 0), the discs are both mas- 
sive and at a high mass-fraction with respect to the star, so 
growth into larger bodies could even start at this very early 
phase when the star is only part-completed. In conclusion, 
most of the dust discs seen around Class II classical 
T Tauri stars are too insubstantial for commencing 
to grow planets, and so these should be relics, in the 
sense that the initial mass reservoir in dust grains 
has already been converted to larger forms with lit- 
tle millimetre emission. This growth process around 
a particular star could have reached anywhere from 
small rocks up to planetesimals or even completed 
gas giant planets, especially if the rapid disc insta- 
bility mechanism operates. 
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